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Abstract 

The Pauli matrices can be derived from the Boltzmann treatment of a two 
dimensional electron gas generated from two different crystal structures. The 
locally created electron gas in A and B wing of a Bell-type experiment enables 
entangled classical currents. We argue for a micro Boltzmann distribution 
■ interpretation of the EPR paradox. 
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1. INTRODUCTION 



Providing a physical theory for local hidden variables [l| is difficult. If, 
however, a macroscopic (semi) classical treatment arrives at entanglement, 
this treatment can model local hidden variables in the micro domain. In 
this paper we will focus the attention on a two-dimensional electron gas. 
The distribution of a two dimensional electron gas can be obtained from the 
semiclassical Boltzmann electron transport theory. 

1.1. Boltzmann transport 

The Fermi-Dirac distribution of electrons in a two dimensional electron 
gas without electric and magnetic fields is ji[ 

/o(k) = ] r ^— t (l) 



1 + exp 



k B T 



Here, k = (k x ,k y ) is the wave vector, k^ is the Boltzmann constant, T the 
absolute temperature and /i the chemical potential. If an electromagnetic 
source is available the acceleration ^v(k) of the electrons is influenced by 
the Lorentz force. 

iv(k) = ~(E(k)+v(k)xB(k)) (2) 
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The acceleration is related to velocity in the wave vector -Ar^rr = -4rk. m* 

J m* at m* 

is the effective mass of the electron. 
1.2. Current 

The distribution / = /(k, t) now differs from equation (TjQ). It is still 
assumed that /(k, t) = /(k + kd£, £ + <it). If the difference between /(k, t) 
and /o(k) is denoted by g(k,t) the current can be approximated by 

j(*) = -2 y 7S2 ev(k)^(k, (3) 

The 2 in equation fl3]) arises from electron spin (considered phenomenolog- 
ically). Note k ^ 0. The linear Boltzmann equation for g(k, t) is equal 
to 

^| = (V k / (k)) ■ |E(k) + (V k g(k, t)) • £v(k) x B(k) (4) 

In a zero magnetic field or when Vk<?(k, t) ~ v(k) + B(k) the linear Boltz- 
mann equation can be solved by 

g(k,t) = Ar(k)|(V k / (k))-E(k) (5) 



Ar(k) = r(k) — to the relative relaxation time. In the general theory |2| 
V k /o = %§fev(k) with 

lim = -5 (£ - Ep) (6) 



Because of non-zero velocities, fl6]) holds approximately [2|. The Ep is the 
Fermi energy. Given, j = crE the conductance 2x2 matrix can be written 

° = 2^ / d2k6 (E(k) ~~ Ef) Ar ( k ) v ( k ) ® v ( k ) ( ? ) 

The tensor product <8> is: (x £g> y)jj = for z, j = 1,2 and x = (a?i, X2) 
similarly y = (yi,y 2 ). 
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2. PAULI MATRICES 

Subsequently, the wave vector k is transformed with k x = kcos((p) and 
k y = ksm(<p). The Jacobian is k = yjk'* + k*. Hence 

e 2 r 

a = — — / kdk dtp 5 (E(k, <p) - E F ) Ar(k, ip)v(k, tp) <g> v(k, tp) (8) 
2nir 2 J 

Note that, because k ^ 0, we also see k ^ and/or dp ^ 0. The next step is 
to observe that E(k,cp) in (jSJ) depends on the velocity v(k,tp) (see 0]). The 
subsequent transformation is to write = E(k,cp) = |m*v 2 (/c, <£>) together 
with (ft = (p. For completeness, = fc(i£, 0) and y> = tp(E, <p). We do not use 
quadratic dispersion for E. The Jacobian for the transformation is 

dk/dE dk/d(f) 
dtp/dE d<p/d(j) 

With this particular transformation dip/dE = and dip /dip = 1. Let us take 
radial kinetic energy only E = |m*v 2 (fc). Hence, 

With ||.|| the Euclidean norm like e.g. in ||k|| = k = v/fc|~+A^. Hence, 
the Jacobian in (j9]) is J = ty(.E') = dk/dE which, in k and looks like 
MlvWII^jp}" 1 - Generally, the ||v(ife,^)|| > 0, despite T -> 0. The 
conductance in (jSj) now can be transformed into 

2 rOO /•+7T 

a = dE d(j)S(E - E F )Ar'(E,(f))v(E,<J))®v(E,(J)) (11) 



J 



(9) 



with Ar f (E,(f)) = k(E,<p)Ar(E,(j))w(E). Subsequent integration over E 
(Ep > 0) results into 

The v generally depend on E and <fi. The first step to the Pauli matrices for 
(IT2"|) is to transform the velocity vector to an 'associated' form: u(Ep, 0) = 
v( J E F ,0) v /Ar / ( J B F ,0). Note that even when v is independent of the u can 
vary with because of 4> dependece in the relaxation time, u can be real or 
imaginary. We have Ar'(E,(f)) = k(E , <p)Ar(E , <f>)w(E) , with real positive 
or negative, relative relaxation time At(E, 0) and, real positive or negative, 
Jacobian w(E). Let us specify Ui(E F , <ft) for i — 1, 2 and derive the first Pauli 
matrix. 
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2.1. Pauli's a x for conductance 

Suppose, the Heaviside H is denned by, H(x) = 1 for all x > and 
H(x) = for all x < 0. Moreover, suppose a G (0, f). Then for G (— a, a), 



8u x> i = Ui(E F , 0) = k(E f )^^ v / tan(0)if(0 + a)H(a - 0) 
Su x ,2 = u 2 (E F ,(j)) = k(E f )J ±^Jcot(4>)H((j) + a)H(a - 0) 

Note that in this definition, suppressing the E F dependence notation for 
the moment, we see uf(0) < when < and u 2 (<p) > when > 
0. The indication Su Xt i refers to the choice of relatively small changes in 
associated velocity for the Pauli matrix a x . For a given > and G 
(0, a) we assume in approximation the conservation of total kinetic energy 
for 'differential associated velocities' 5u x (E,(p) = u(E,<f)). I.e. the changes 
in kinetic energies when electron and hole are created occur 'balanced' in the 
u. 

l -m* 1 1 u(E, 0) 1 1 2 + l -m* \\\x{E, -0) 1 1 2 = (14) 

It is believed that electron escape from a crystal structure leaving behind a 
hole can be pictured in equation (|T5|) . Note that equation f[T4"j) is valid in u 
not in v terms. Integrating for u\ for instance, using (I12p then (punching 
= 0) gives 



e 2 



Or, equally 



= ^^^ / + "^ tan ^) = - C [ lo S I cos (^)l]-« ( 16 ) 

with C = ^§^2^2- From ([16]) follows ai 5 i = — C{log | cos(a)|— log | cos(— a)|}. 
Hence, cr^i = 0. A similar argument applies to 0*2,2 where the integration 
shows that <t 2)2 = — C{log | sin(a)| — log | sin(— a)\} = 0. If we subsequently 
turn to o~\ 2 = o~ 2 1 then we see 



Cl,2 



k 2 (E f ) e 2 f +a J± , —r , -r e 2 K 2 (E F ) 



4 



Hence, with n(Ep) = ZLy ^ the Pauli matrix a x can be obtained for the 
conductance matrix. However, because of infinitesimal changes in 'associated 
velocity' (Tl3|) . the factor is maintained. Hence, 

_ _ e 2 n\E F ) ( 1 

2->x - 



2hn 2 \ 1 

k sufficiently small. Note that 'going through a singularity' for the cot inte- 
gral provides a zero result of the integral when the = is 'cut out' of the 
integration by left-hand — e and right-hand e for e — > + . This is so because 
| sin(— e)| equals | sin(e)|. Hence, the integral of both the tan as well as the 
cot containing associated 'velocities' in equation ( TT3l) can be given by the 
e — > + sum of f_ e and f n integration operations and is written *". The 

Heavisides in (fl3|) convert the integrals into f* a . 

2.2. (7 y conductance 

The o y Pauli matrix can be similarly derived for the conductance. Let us 
define the 'associated velocity' entries 

Suy A = Ul (E F , 0) = X(Ep) cos(0/2)#(0 - P)H(\ - 0) , , 

5u y , 2 = u 2 {E F , 0) = % \{Ep) sin(0/2)#(0 - f )H (tt - - 0) {lJ) 

with f3 G (0, a). In the first place we may note that, dissimilar to the a x case 
in secion - I2.1[ the canceling of 'associated' kinetic energy terms occur in 
a non-symmetrical way: i.e. |m*||u(0)|| 2 + |m*||u(0 / )|| 2 = where (0',0) 
solves cos 2 (0) -sin 2 (0') = with G (/3, f ) and 0' G (|,7r — /3). If the aim is 
to explain the geometric asymmetry in electron - hole creation geometry then 
perhaps the crystal structure, i.e. phonon hindrance can be held account- 
able for the asymmetry. Subsequently, the integral for ui(Ef,4>)u2(E f ,(P) 
vanishes because the domains are disjoint. For the er^i term we see 

0-1,1 = 2 ^ 2 F) J cos 2 (0/2)# (20) 

Because cos 2 (0/2) = |(1 + cos(0)) it follows that 

°~i,i = 6 (f — $ + 1 ~~ srn (/3)}- Similarly for cr 2 ,2 we may derive 
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Because sin 2 (0/2) = ~(1 — cos(0)) it follows that 

cr 2 ,2 = — e {f ~ + 1 ~ sin(/3)|. Hence, when we take X(E F ) and 

k(E f ) such that 

X 2 (E F ) rvr 



then we find 



- 2 (^) = ^P{^-/3 + l-sin(/3)} (22) 



_ e V(iS F ) MO 

3. ENTANGLED CURRENTS 

Suppose in experiment, similar to Bell's jsd, jH we have a distant A and 
_B wing where, at both sides, a localized but identical electric field with vecor 

7^ 2/wr 2 * T 2H-K 2 , _ _ . . . 

E = 2 2 / F X E = 2 2^ \ ( 1 ? E2) 24 

e z K z {E F ) e z K z (E F ) 

reigns. Furthermore we insist that E 2 = E\ + E\ = 1. If we take the two 
'crystal structures' we talked about in the previous two sections and generate 
electron-hole pairs, then at first instance, electrons in the 2D gas from the 
horizontal H x related crystal and from the vertical S y related crystal co-occur 
in the anglular interval (/3,a). In this area per wing one can 'mix the two 
conductances' and £ y with a unit parameter vector a = (01,02) and/or 
b = (61, 6 2 )- E.g. Ea(oi, a 2 ) = a{E x + a 2 Ti y . The mixed current at A is then 
equal to = S j4 (ai,a 2 )E. Similarly, one can define a T, B (bi,b 2 ) such that 
Jb = £b(£>i, b 2 )E. Note that we can take a 2 x 100% to indicate the percentage 
at the A wing of S^-crystal electrons and similarly for the a 2 , x 100%. The 
same thing can be supposed for the entries of the b vector. Furthermore, the 
mixing percentages are transformed into angles. We suppose, a\ = cos(iPa) 
and d2 = sin (^4) together with b\ = cos(iJjb), b 2 = sin^s) and project the 
if) in an interval [7, At] C (0, |). Subsequently, we may note that because of 
fl24|) and the expressions for S x and T> y in resp. ( ITS]) and fl23|) the following 
inner-product for and js employs Pauli matrices as in a quantum corelation 
(see e.g. Q). 

j B = E T [cos^aVz + sin(^ A )cr y ] [cos(^b)ct x + sin(V>B)er v ] E (25) 
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Note that, a J = a x and a J = a y together with o 2 x = <jy = l 2 x2- If we in 
addition inspect e.g. the term cos(ipA) sin(^B)E T cr x cr y E it is noted that this 

( ~E 2 \ 

equals cos(^a) shifts) E 2 ) ( ^ J = Hence, the inner product of 

the A and B wing current is js = cos(ipA — iPb)- The inner product of the 
two current vectors are entangled. 

4. CONCLUSION 

It was demonstrated that under classical local conditions, entangled mea- 
surements can be simulated that commonly are thought only to occur in 
non-local quantum theory. Let us summarize the experiment: The source 
in a Bell-type experiment can be the cloning and subsequent sending in two 
different directions of the electric field vector E. Two pairs of crystals are 
employed to generate the mixing of conductances S x and S y , in the A and 
conductances H x and £ y , in the B wing. The electrons escape from the 
crystal surfaces (E x parallel x + axis and parallel y + axis and xl.y) and 
locally mix in G (/?,«). If the current vectors that are created from the 
'cloned' electric fields vectors are transported to a measuring system O a 
current-current entangled inner product ]b = a • b can be observed from 

= [a x a x + a 2 a y ]E A and j B = [hcr x + b 2 cr y ]E B . Hence, j B can in princi- 
ple violate the CHSH but note that the current vectors are created by local 
means. The parameters a and b refer to mixing percentages of electrons with 
<p G a): i.e. in the A wing a\ x 100% of the electrons from the E x crystal 
and a\ x 100% from the S y crystal and similarly in the B wing. Perhaps 
one would like to argue against entanglement in currents but the outcome in 
the innerproduct j# is the same as quantum mechanically. Entanglement 
refers to something unobservable and obtains its meaning from its use [5] . Es- 
sentially it is concluded from the experimental correlation. Unless we refered 
unwittingly to macroscopic quantum electron gas, Bell's correlation needs 
not be nonlocal in its origin. It can be claimed that the structure laid down 
here is the physics of locality and causality that Einstein had in mind. The 
question then of course, similar to Madelung's hydrodynamic interpreation 
of qm 0], is: 'if referred to the microphysics domain, what does the distri- 
bution in the Boltzmann equation distribute'. On the mathematical side 
the result is in accordance with 0], [ij], 10 . 



An additional point arises if we note that mathematically there is a con- 
finement to a mixing interval (f> G (ft, a), with > 0. It can be that, macro- 
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scopically, the mathematical structure does not violate the CHSH because 
of the mixing angle restriction. This does not disqualify the mathematical 
structure used for local hidden variable purposes but shows the difference of 
macro vs micro Boltzmann distributions. 
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